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We analyze photo-refractive incoherent soliton beams and their interactions in Kerr-like nonlinear 
media. The field in each of Ad incoherently interacting components is calculated using an integrable 
set of coupled nonlinear Schrodinger equations. In particular, we obtain a general A^-soliton solution, 
describing propagation of multi-soliton complexes and their collisions. The analysis shows that the 
evolution of such higher-order soliton beams is determined by coherent and incoherent contributions 
from fundamental solitons. Common features and differences between these internal interactions are 
revealed and illustrated by numerical examples. 



One of the most noted discoveries of modern soliton 
science is that solitons can be excited by an incandes- 
cent light bulb instead of a high power laser source [Q . 
This produces "incoherent solitons" they can exist 

in photorefractive materials which require amazingly low 
powers to observe highly nonlinear phenomena It 
is also remarkable that, in certain conditions, incoherent 
solitons in photorefractive materials can be studied using 
coupled nonlinear Schrodinger equations (NLSE) . 

In general, coupled nonlinear Schrodinger equations 
(NLSE) can be applied to various phenomena. These 
include incoherent solitons in photo-refractive materi- 
als, plasma waves in random phase approximation 
multicomponent Bose-Einstein condensate ||] and self- 
confinement of multimode optical pulses in a glass fiber 
10 . Therefore its solutions are of great interest for the- 
oretical physicists. In special cases these equations are 
found to be integrable ijlll . Then, in analogy with single 
(scalar) NLSE (when the number of equations, M, is 1) 
^ and the Manakov case {M = 2), the total solu- 
tion consists of a finite number (N) of solitons and small 
amplitude radiation waves. The former is defined by the 
discrete spectrum of linear {L,A) operators [p"2|jl^ ] and 
the latter is defined by the continuous spectrum. Most 
applications deal with the soliton part of the solution as 
it contains the most important features of the problem. 
Moreover, a localized superposition of fundamental soli- 
tons can be called "multisoliton complex" . An incoherent 
soliton is a particular example of a multisoliton complex 

The cases M — 1 and M = 2 have been extensively 
discussed in the literature [ p2|jr^ ]. On the other hand, 
results for general M are scarce. The linear (L, A) op- 
erators are important elements for the inverse scattering 
technique, which can be considered as a basis for integra- 
bility of M coupled NLSEs. Moreover, it has been shown 
| p5[ that A^-soliton solutions of M coupled NLSE can 
be found using a simple technique which is an extension 
of the theory of reflectionless potentials In recent 



works ||lj,|l3l cases when each component has only one 
fundamental soliton have been considered. It was demon- 
strated that, in this configuration, the formation of sta- 
tionary complexes may be observed, and corresponding 
solutions for M = iV < 4 were presented in explicit form 

So far, only the case of complete mutual incoherence 
of the fundamental solitons has been considered. In this 
case the multisoliton complex can also be viewed as a 
self- induced multimode waveguide ||l^. The general case, 
where fundamental solitons in the multisoliton complex 
interact both coherently and incoherently, has not been 
analyzed. Such interactions may be observed if N is 
larger than M, so that each component has not less than 
one fundamental soliton. In general, each fundamental 
soliton can be "spread out" among several components. 
We will refer to this effect as mixed "polarization" of 
fundamental solitons. However, in order to capture dis- 
tinctive features of coherent and incoherent soliton inter- 
actions, we will focus on a special case which is important 
for incoherent solitons. Specifically, we consider a situa- 
tion where all the fundamental soliton polarizations are 
mutually parallel or orthogonal, and thus are conserved 
in collisions jlj] . Due to the symmetry of the NLSE with 
respect to rotations in functional space, hereafter we as- 
sume that each fundamental soliton is polarized in one 
component only. It is for this case that we present new 
explicit iV-soliton solutions of M coupled NLSEs, and we 
discuss the new physics which it brings into the theory. 

We consider propagation of an incoherent self-trapped 
beam in a slow Kerr-like medium and write the set of 
coupled NLSEs in the form WS^: 
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where ipm denotes the m-th component of the beam, z 
is the coordinate along the direction of propagation, x is 
the transverse coordinate, and 
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is the change in refractive index profile created by all in- 
coherent components of the light beam, where the am 
(> 0) are the coefficients representing the strength of the 
nonlinearity, and M is the number of components. 

Solutions in the form of multisoliton complexes of 
Eq. (|^) and their collisions can be obtained using the 
formalism of |l5| , [l^ with some refinements. First, we 
introduce functions Uj{x, z) as solutions of the following 
set of equations: 
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where A'' is a total number of fundamental solitons, 
Cj = Xj exp {kjXj + ikjZj/2), xj — x — xj and zj — z — Zj 
are shifted coordinates, and Xj ^re arbitrary coefficients. 
The values Xj and Zj characterize the initial positions of 
fundamental solitons, but the actual beam trajectories 
may not follow the specified points due to mutual inter- 
actions between fundamental solitons. Each fundamental 
soliton is characterized by an eigenvalue kj — Tj ~\- i^j. 
Its real part, rj, determines the amplitude of the fun- 
damental soliton, while the imaginary part, iij — tanOj, 
accounts for the soliton velocity (i.e. motion in trans- 
verse direction) . Here 9j is the angle of the fundamental 
soliton propagation relative to the z axis. 

To distinguish coherent and incoherent contributions 
to the multi-soliton complex, we use variables nj, which 
represent the number of the component where the j-th 
soliton is located. Thus, two fundamental solitons with 
rij = Tim are coherent, and they are incoherent otherwise. 
Now we can write the expression for the matrix D: 
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Finally, the A^-soliton solution of the original Eq. (|^) 
can be obtained by adding up of all the Uj corresponding 
to a given component number m: 
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Note that the number of terms in the sum is exactly the 
number of fundamental solitons polarized in this compo- 
nent, viz. Njn, and the total N is X]m=i -^m- 

One of the features of this approach is that coherent 
fundamental solitons are "split" among all the uj func- 
tions for a given component. However, when obtaining 
analytical solutions in explicit form, it is possible to sepa- 
rate fundamental solitons by combining terms with corre- 
sponding propagation constants. Consequently, we write 
the exact solutions for a different set of functions uj , with 



each of them containing one fundamental soliton (at dis- 
tances where coherent interactions are small) . These are 
combined into the original functions in the following way: 



I J J . I L -I lit J ' 

The coefficients Xj are arbitrary, and we can choose 
particular values for them: 
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where bjm = {kj + k'^)/{kj — km), and the square root 
value is taken on the branch with positive real part. This 
step significantly simplifies further analysis, as the result- 
ing solution will acquire a highly symmetric form. 

Finally, the explicit expressions for solutions can be 
found as sums over specific permutations: 
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Here L denotes four sets of indices (Li, L2, ^3, L4). The 
summation is performed over all combinations in which 
the given set of soliton numbers (for example, {1, . . . , N}) 
can be split among all the Lj . When performing permu- 
tations, Li, L2 are only filled with numbers of mutually 
coherent solitons (thus the number of elements in these 
sets is the same). 

The coefficients and functions from (|^) are determined 
for each realization of the permutation L as follows: 
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Fl{x,z)= cos(S'g) cos(S't) cosh(S'b) — 
sin(S'g) sin(S'f) sinh(S'b), 

— I -L; -'c-'sg -'sg -'mg-tg-'sb -'sb -'mb-'mgb, 

Fl{x,z)= cos(S'|)cos(S'f)cosh(5;^)- 
sin(S'l) sin(5f) sinh(S'j^). 

Here we used \Li \ to denote the number of elements in the 
set. Note that the F functions are written in the simplest 
form in terms of trigonometric and hyperbolic functions, 
due to the specific choice of coefficients in Eq. (H) . 

The variables introduced above are the following sums 
and products over the Lj sets: 
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Here the "U" operator is used to merge the sets, 
and the variables + ijj — kjXj + ik'^Zj/2 (with 



(3j and jj real), rjj, 



logi\ik,-k,^){k,+k*J\)/2, 



sign I TT 



= \bjm\, Vjni = arg(l/6j„)/2, Sjr 



arg 



)]}sign (m-j), 

&j7n = (fcj +k^)/{kj — km)- The function arg is supposed 
to give values in the interval [0, 27r), and 



sign 



1, a;>0, 
-1, a; < 0. 



Note that only [3j and 7-,- depend on the coordinates 
(x, z). All the other coefhcients are expressed in terms 
of the wave numbers kj and constant shifts in positions 
(xj, Zj) of the TV fundamental solitons. As the total so- 
lution has translational symmetry, one of the shifts can 
be fixed, so that the number of independent parameters 
controlling the multisoliton complex is 2N — 1. 

If an incoherent soliton consists only of orthogonally 
polarized fundamental solitons {uj = j, N = M), and all 
are propagating in the same direction, then its transverse 
intensity profile remains stationary Jl^ . In this partic- 
ular case, the general expressions (|8|)are radically sim- 
plified, since, due to the above-mentioned restrictions on 
the permutations, the sets Li and L2 are always empty. 
Hence, we obtain: 
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Fl = cosh(5b), fI = coshiSi). 

Note that here we have neglected a common multiplier in 
Cl and C;^, as these coefficients determine respectively 
the denominator and numerator in the expression for Uj . 




FIG. 1. Stationary propagation of an incoherent soliton 
consisting of eight completely incoherent fundamental solitons 
(polarized in different components). 

Now we present numerical examples to illustrate these 
results. An example of a stationary incoherent soliton 
consisting of eight components {N — Ad — 8) is shown 
in Fig. 1^. The profiles of the constituent fundamental 
solitons, and their superposition as a whole, are deter- 
mined by the wave numbers and relative shifts along the 
X axis. In this configuration, the shifts in propagation di- 
rection, Zj, correspond to arbitrary phase changes of dif- 
ferent components, but these do not influence the evolu- 
tion due to the incoherent nature of the inter-component 
interactions. 




FIG. 2. Evolution of an incoherent soliton with multi-scale 
periodic "beating" due to internal coherent interactions (8 
fundamental solitons in 5 components). 

On the other hand, if > M, two or more of the fun- 
damental solitons are polarized in the same components, 
and thus interact coherently. If the inclination angles of 
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the fundamental solitons are all the same, the beam will 
remain localized upon propagation. Such a multi-soliton 
complex is an incoherent soliton with an intensity profile 
which evolves periodically or quasi-periodically, as shown 
in Fig. |2[ These oscillations, appearing due to internal 
coherent intra-component interactions, are a general fea- 
ture of incoherent solitons, and can be eliminated only 
in specific cases, as discussed earlier. It follows that spa- 
tial "beating" always accompanies the interaction of fun- 
damental solitons of a single NLSE, which agrees with 
previous studies . 

Our explicit solution (|^) also describes collisions of in- 
coherent solitons. As mentioned earlier, the polarizations 
of the fundamental solitons are preserved in collisions 
(provided they are orthogonal or parallel), and thus the 
degree of internal coherence doesn't change. However, 
the shifts of the fundamental soliton trajectories differ, 
and this results in the incoherent solitons changing their 
shapes. These transformations can be seen clearly in 
Fig. I 




FIG. 3. Collision of a completely incoherent soliton (con- 
sisting of two orthogonally polarized fundamental solitons) 
and an incoherent soliton with internal coherent contributions 
(6 fundamental solitons in 5 components). 

To calculate the shifts, we use the fact that in the 
expression for soliton profiles, Uj, given by Eq. (^, the 
denominator U is real, and the numerator does not de- 
pend on the coordinates of the corresponding fundamen- 
tal soliton {xj,Zj). It is then straightforward to take 
appropriate limits and calculate the shift of j-th funda- 
mental soliton along the x axis due to collisions: 

^ — S^ ±1 21n(cj„i), nj = n„. 

Here the summation involves the fundamental solitons 
which feature in the collisions. The "-I-" sign corresponds 
to the case when colliding soliton number m comes from 
the right (i.e. has larger x coordinate before the impact), 
and the "— " sign when from the left. This is a general- 
ization of the expressions found in llj] . 



In summary, we have obtained a general iV-soliton 
solution of M coupled nonlinear Schrodinger equations 
which describes multi-soliton complexes supported by 
a Kerr-type nonlinearity. A particular example is an 
incoherent soliton in a photo-refractive medium. We 
have revealed that the properties of multi-soliton com- 
plexes, which are superpositions of fundamental soli- 
tons with orthogonal or parallel polarizations, are de- 
termined by internal interactions, both phase-insensitive 
inter-component and coherent intra-component, with the 
latter resulting in spatial " beating" . Using our exact re- 
sult, we also analyzed collisions of incoherent solitons. 
We showed that the re-shaping of incoherent solitons af- 
ter collisions are characterized by the relative shifts of 
the fundamental solitons, and these are calculated using 
a simple analytical formula. These distinctive features of 
incoherent solitons are illustrated by numerical examples. 

The authors are part of the Australian Photonics CRC. 
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this manuscript. 



[1] M. Mitchell and M. Segev, Nature (London) 387, 880 
(1997). 

[2] M. Mitchell, Z. Chen, M. Shih, and M. Segev, Phys. Rev. 

Lett. 77, 490 (1996). 
[3] M. I. Carvalho, T. H. Coskun, D. N. Christodoulides, M. 

Mitchell, and M. Segev, Phys. Rev. E 59, 1193 (1999). 
[4] G. Duree et al, Phys. Rev. Lett. 71, 533 (1993); 
[5] A. A. Zozulya, D. Z. Anderson, A. V. Mamaev, and M. 

Saffman, Europhys. Lett. 36, 419 (1996). 
[6] M. D. Iturbe-Castillo et al, Appl. Phys. Lett. 64, 408 

(1994) 

[7] V. Kutuzov, V. M. Petnikova, V. V. Shuvalov, and V. A. 

Vysloukh, Phys. Rev. E 57, 6056 (1998). 
[8] A. Hasegawa, Physics of Fluids 20, 2155 (1977). 
[9] E. P. Bashkin and A. V. Vagov, Phys. Rev. B 56, 6207 

(1997). 

[10] A. Hasegawa, Opt. Lett. 5, 416 (1980). 

[11] V. G. Makhan'kov and O. K. Pashaev, Teor. Mat. Fiz. 

53, 55 (1982) [Theor. Math. Phys. 53, 979 (1982)]. 
[12] V. E. Zakharov and A. B. Shabat, Zh. Eksp. Teor. Fiz. 

61, 118 (1971) [Sov. Phys. JETP 34, 62 (1971)]. 
[13] C. B. Manakov, Zh. Teor. Eksp. Fiz. 65, 505 (1973) [Sov. 

Phys. JETP 38, 248 (1974)]. 
[14] N. N. Akhmediev, W. Krolikowski, and A. W. Snyder, 

Phys. Rev. Lett. 81, 4632 (1998). 
[15] Y. Nogami and C. S. Warke, Phys. Lett. 59A, 251 (1976). 
[16] I. Kay and H. E. Moses, J. Appl. Phys. 27, 1503 (1956). 
[17] A. Ankiewicz, W. Krolikowski, and N. N. Akhmediev, 

Phys. Rev. E (to be published in May, 1999). 
[18] C. S. Gardner, J. M. Greene, M. D. Kruskal, and R. M. 

Miura, Comm. Pure Appl. Math. 27, 97 (1974). 
[19] N. Akhmediev and A. Ankiewicz, Solitons, nonlinear 

pulses and beams (Chapman & Hall, 1997). 



4 



